Large Fluctuations in Driven Dissipative Media 
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We analyze the fluctuations of the dissipated energy in a simple and general model where dissipa- 
tion, diffusion and driving are the key ingredients. The full dissipation distribution, which follows 
from hydrodynamic fluctuation theory, shows non-Gaussian tails and no negative branch, thus vi- 
olating the fluctuation theorem as expected from the irreversibility of the dynamics. It exhibits 
simple scaling forms in the weak- and strong-dissipation limits, with large fluctuations favored in 
the former case but strongly suppressed in the latter. The typical path associated to a given dissipa- 
tion fluctuation is also analyzed in detail. Our results, confirmed in extensive simulations, strongly 
support the validity of hydrodynamic fluctuation theory to describe fluctuating behavior in driven 
dissipative media. 



Fluctuations are inherent in most physical processes, 
hinting at the hectic microscopic dynamics that results 
in observed macroscopic behavior. They encode essen- 
tial physical information despite their apparent random 
origin, and their investigation is opening new paths to un- 
derstand physics far from equilibrium [TH5]. The study 
of fluctuation statistics of macroscopic observables pro- 
vides an alternative way to derive thermodynamic poten- 
tials, a complementary approach to the usual ensemble 
description. This observation, valid both in equilibrium 
and nonequilibrium [3J [3], is most relevant in the lat- 
ter case where no bottom-up approach exists yet con- 
necting microscopic dynamics with macroscopic proper- 
ties. The large deviation function (LDF) controlling the 
statistics of these fluctuations [1] plays in nonequilibrium 
physics a role akin to the equilibrium free energy. Key 
for this emerging paradigm is the identification of the 
relevant macroscopic observables characterizing nonequi- 
librium behavior. The system of interest often conserves 
locally some magnitude (a density of particles, energy, 
momentum, charge, etc.), and the essential nonequilib- 
rium observable is hence the current or flux the system 
sustains when subject to, e.g., boundary-induced gradi- 
ents or external fields. In this way, the understanding 
of current statistics in terms of microscopic dynamics 
has become one of the main objectives of nonequilibrium 
statistical physics, triggering an enormous research effort 
which has led to some remarkable results jlHH] . These re- 
cent advances are however restricted to nonequilibrium 
conservative systems characterized by currents. On the 
other hand, many nonequilibrium systems are character- 
ized by an irreversible dissipation of energy and need a 
continuous input of energy in order to stay stationary. 
Thus the relevant macroscopic observable here is not only 
the current but also the dissipated energy. This broad 
class of systems includes granular media, chemical reac- 
tion systems, dissipative electronic and biophysical sys- 
tems, turbulent fiuids, and in general all sort of reaction- 
diffusion systems where dissipation, diffusion and driving 
are the main ingredients. Fluctuations in dissipative me- 



dia have been much less studied, most probably because 
its physics is far more complicated as a result of the ir- 
reversibility of their microscopic dynamics. 

In this paper we address this issue by studying fiuctu- 
ations of the dissipated energy in a paradigmatic model 
of driven dissipative media, using both hydrodynamic 
fluctuation theory (HFT) and extensive numerical sim- 
ulations. We consider a generalization of the Kipnis- 
Marchioro-Presutti (KMP) model of heat transport [5] 
to dissipative media. The system is deflned on a one- 
dimensional (ID) lattice of size N, where each node 
i e [1, A^] is characterized by an energy pi. Dynamics is 
stochastic and proceeds via collisions between randomly- 
chosen nearest neighbors + 1) such that a fraction 
1 — a of the total pair energy = + p^+i is dissi- 
pated out of the system, and the remaining energy aYn is 
randomly redistributed within the pair. Thus, the post- 
collision energies are p'^ = zaYn and p[j^i = (1 — z)aY,i, 
with z a random number uniformly distributed in [0, 1]. 
In addition, boundary sites (i = 1,A'^) may interchange 
energy with a heat bath at temperature T so the sys- 
tem reaches a steady state where energy injection and 
dissipation balance each other. For conservative dynam- 
ics (a = 1), this model represents at a coarse-grained 
level the physics of many diffusive systems of technolog- 
ical as well as theoretical interest, and plays a main role 
in nonequilibrium statistical physics as a benchmark to 
test theoretical advances [IHS]. Our generalization, see 
also Hn] , contains the essential ingredients characterizing 
most dissipative media, namely: (i) diffusive dynamics, 
(ii) bulk dissipation, and (iii) boundary injection. More- 
over, our system can be regarded as a toy model for a 
dense ID granular gas, in which particles cannot freely 
move but can collide with their nearest neighbours, los- 
ing a fraction of the pair energy and exchanging the rest 
thereof randomly. The parameter a can be thus con- 
sidered as the analogue to the restitution coefficient in 
granular systems This model is an optimal candi- 

date to study dissipation statistics because: (a) one can 
obtain explicit predictions for the LDF, and (b) its sim- 
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pie dynamical rules allow for a detailed numerical study. 
The chances are that our results remain valid for more 
complex dissipative media described at the mesoscopic 
level by a similar evolution equation. We report below 
simulation results for the statistics of the dissipated en- 
ergy in this model using both standard simulations and 
an advanced Monte Carlo method [12] . The latter allows 
the sampling of the tails of the distribution, and implies 
simulating a large number M of clones of the system. 

Our model belongs to a large class of dissipative media 
characterized at the mesoscale by the following rescaled 
Langevin equation (T^] 

dtp{x, t) = -dxj{x, t) - iyp{x, t) . (1) 

Here t) is the energy density field, with x e 
[— |, \] and boundary conditions p(±^,i) — T, j{x,t) = 
Q[p{x,t)] + £.{x,t) is the fluctuating current, with lo- 
cal average behavior given by Fourier's law, Q[p\ = 
—D[p]dxp{x, t), D[p\ is the diffusivity functional, and v is 
the mesoscopic dissipation cocfhcicnt. The first term on 
the rhs of eq. ([T]) describes the diffusive spreading of en- 
ergy, while the second term gives the rate of energy dissi- 
pation in the bulk. The current noise term ^{x, t), which 
accounts for microscopic random fluctuations at the 
mesoscopic level, is Gaussian and white with {^{x, t)) =^ 
and {C{x,t)^{x',t')) = N-^a[p]S{x - x')S{t - t'), be- 
ing a[p] the mobility functional. This Gaussian fluctu- 
ating field is expected to emerge for most situations in 
the appropriate mesoscopic limit as a result of a central 
limit theorem: although microscopic interactions can be 
highly complicated, the ensuing fluctuations of the slow 
hydrodynamic fields result from the sum of an enormous 
amount of random events at the microscale which give 
rise to Gaussian statistics of 0{N^^/^) at the mesoscale. 

The probability of observing a particular history 
{p{x,t), j{x,t)}Q of duration r for the density and cur- 
rent fields follows from a path integral and obeys a large 
deviation principle Ij, P{{p,j}o) ^ exp{+NXr[p,j]}, 
with a rate functional 
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with p{x,t) and j{x,t) coupled via the balance equa- 
tion ([!]). Eq. ([2| expresses the Gaussian nature of lo- 
cal current fluctuations around its average (Fourier's) 
behavior, and quite remarkably it takes the same form 
as in many conservative diffusive systems This 
stems from the subdominant role of the noise affecting 
the dissipative term in eq. which scales as N~'^/'^ 
and thus is negligible against the current noise in the 
mesoscopic limit J13j. We are interested in the prob- 
ability of a fluctuation of the time-averaged dissipated 
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FIG. 1: (Color online) Scaling of the dissipation LDF in the 
quasi-elastic limit {f ^ 1) for A'^ = 50, T = 1 and varying v. 
The solid and dashed lines are the HFT prediction and Gaus- 
sian estimation, respectively. Small points around the peak 
were obtained in standard simulations. Inset: Gonvergence of 
time-averaged dissipation and sketch of the probability con- 
centration associated to the large deviation principle. 



[T] and the dissipation LDF G{d) is related to Xr[p,j] 
via a saddle-point calculation in the long time limit, 
G{d) — maXpjX7-[p, j], such that the optimal profiles 
Pq{x, t; d) and jo{x, t; d) solution of this variational prob- 
lem are consistent with the value of d and the balance 
equation ([T]). These optimal profiles can be interpreted 
as the ones the system adopts to facilitate a given dissi- 
pation fluctuation. For conservative systems it has been 
shown that optimal profiles equivalent to ours are in fact 
time-independent in a broad regime as a result of an ad- 
ditivity property of fluctuations [11[S1IH], a scenario which 
may break down only for extreme fluctuations |14j . As- 
suming now time-independent profiles in our dissipative 
case, the balance equation simplifies to j' (x) + h'p{x) — 0, 
with the prime denoting spatial derivative, and 



.1/2 li.jix)-D[-^-]f{x) 
G{d) = - min / dx^ . (3) 

The functional above is a generalized Lagrangian and the 
associated optimal current profile jQ(x;d) thus obeys a 
fourth-order Euler-Lagrange equation |13j . with bound- 
ary conditions j'Q{±^;d) = —vT and jo{—^;d) — 
— jo(i?;c?) = d/2. These follow from the relation / = 
— i/p, the definition of the dissipation d and the symme- 
try of our problem around x ~ 0. 

For the KMP model D[p] = 1/2 and a[p] = p"^ 0113]. 
The stationary average solution to eq. ([T]) hence reads 
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FIG. 2: (Color online) Optimal energy profiles for varying 
d/ {(£) measured for u = 10^'' (symbols) and v — 10^^ (dashed 
lines), and HFT predictions (solid Hnes and inset). 
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and the average dissipation is (d) = v J_y^dxpst{x) = 

r-\/2i/ tanh(-y/j//2). Interestingly, these steady profiles 
are similar to the ones found in vibrated granular gases, 
when measured in units of the mean-free-path |15j . The 
dissipation parameter v can be related to KMP micro- 
scopic dynamics via — 2{N + 1)^(1 — a) [13] ■ This 
means that the microscopic dynamics consistent with eq. 
^ must be quasi-elastic (1 — a — 0{i'N~'^)) to guarantee 
that both diffusion and dissipation interplay on the same 
scale at the mesoscopic level [17]. The Euler-Lagrange 
equation for the optimal current profile solution of the 
variational problem eq. ([3| for KMP is 

(5) 



12i^^j(j'^ - 33") + J"(3j" - 4j'j"') + = . 



For small fluctuations of the dissipation, e = (d — 
{d))/{d) <C 1, a regular perturbative analysis of eq. ([s]) 
around the steady-state profile eq. ^ yields G{d) w 
— e^/2A^, i.e. Gaussian statistics for small fluctuations 
around the average dissipation, as expected from the cen- 
tral limit theorem. The variance is = [sinh(2\/2i^) — 
2v^]/[4v^sinh^(V2z^)]. Interestingly, ~ 1/3 inde- 
pendent of V for 1/ <C 1, anticipating the existence of a 
simple scaling form for G{d) in this quasi-elastic limit 
(see below). On the other hand, A^ ~ {2^/2v)^^ for 
V ^ I, suggesting a suppression of dissipation fluctua- 
tions in the strongly inelastic regime. Remarkably, A^ 
turns out to be a good estimate for the variance of the 
dissipated energy, see inset to Fig. [3] 

We now obtain the whole dissipation LDF in the quasi- 
elastic limit V ^ 1. The behavior of the steady-state 
current field and average dissipation for ^ 1, see eq. 
(|4]), suggests the scaling j{x) = wipi^x) and d = v/S. in 
this limit, where il){x) and A (the total energy per site) 
remain of the order of unity. In fact, {d) ~ vT for <C 1, 
i.e. A ~ r. Using this scaling in eq. ([s]) and retaining 
only the lowest order in v we get a differential equation 
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where po(a;) = —'((^^{x), and 6 is a constant given by 

sinh b d 
~ = Jd)- 



(6) 
(7) 



The average behaviour in this quasi-elastic regime is re- 
covered in the limit 6^-0 (i.e. d — (d)). For 6^1 
we have small deviations from the average and Gaussian 
statistics as described above. For arbitrary values of 5, 
eq. ^ gives the optimal profiles for the system to sus- 
tain an arbitrary fluctuation of the dissipated energy d. 
Using these optimal proflles in eq. ([s]) we flnd for z^ ^ 1 



G{d) = 6 tanh 
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with b implicitly given by eq. ([7| in terms of d/{d). 
This means that, in the low dissipation limit ly 1, 
G{d) shows a simple scaling form independent of v when 
plotted against the relative dissipation d/{d), as antic- 
ipated above. The scaling form is fully confirmed in 
Fig. [l] which shows the G{d) as measured for differ- 
ent, small values oi v ^ [10^"^, 10^^] in simulations of the 
dissipative KMP model. The dissipation LDF is highly 
skewed with a sharp decrease for fluctuations d < {d) 
and no negative branch, so fluctuation theorem-type re- 
lations linking the probabilities of a given dissipation d 
and the inverse event —d do not hold [6l [7]. This was 
of course expected from the lack of microreversibility, 
a basic tenet for the fluctuation theorem to apply [16j . 
The limit 6 » 1 corresponds to large dissipation fluc- 
tuations, where G{d) « —^[ln{d/{d))]'^, i.e. a very slow 
decay which shows that such large fluctuations are far 
more probable than expected within Gaussian statistics 
(^ — |(d/(d))2). We also measured the typical energy 
profile associated to a given dissipation fiuctuation, see 
Fig. [2] finding also very good agreement with HFT. Note 
that optimal profiles for varying v 1 also collapse for 
constant d/{d). Moreover, profiles associated to dissi- 
pation fiuctuations above the average exhibit an energy 
overshoot in the bulk. This observation suggests that the 
mechanism responsible for large dissipation fiuctuations 
consists in a continued over-injection of energy from the 
boundary bath, which is transported to and stored in the 
bulk before being dissipated. 

For J/ > 1 we solved numerically eq. (5|) and used 
this optimal profile to obtain G{d). Fig. [3 (top) shows 
the predicted G{d) for increasing, non-perturbative val- 
ues of and results from simulations. A slow but clear 
convergence toward the macroscopic HFT prediction is 
observed as N increases. Such strong finite-size effects 
are expected since the natural lengthscale associated to 
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The reason is that Ref. [T7] studies systems with two 
competing dynamics, one conservative and another non- 
conservative, thus resulting in independent fluctuations 
for the density and dissipation fields. In our theory, as is 
the case in many driven dissipative systems, dissipation 
is linked to the collision process and hence dissipation 
fluctuations are enslaved to density profile deviations. In 
fact, both theories coincide in the limit where the opti- 
mal dissipation profile is given in terms of the optimal 
density field [13]. 

In summary, we have shown that hydrodynamic fluctu- 
ation theory is a powerful tool to describe fluctuating be- 
havior in driven dissipative media. This opens the door 
to further general results in nonequilibrium physics, as 
HFT provides a well-defined scheme based on the knowl- 
edge of the conservation laws governing a system and just 
two transport coefficients. 
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FIG. 3: (Color online) Top: Dissipation LDF measured 
for increasing values of f (sliifted vertically for convenience, 
Criid)) ~ V^). Open symbols correspond to = 50 and 
M = 10^ clones, while for filled symbols iV = 100 and 
M = 10*. Solid lines are HFT predictions obtained by nu- 
merical integration of eq. ([sjl . Notice the (slow) convergence 
toward HFT as A'^ and M increase. Inset: Measured average 
dissipation and its variance as a function of and Gaussian 
estimation. Bottom: Optimal energy profiles measured for 
u = 10, N = 50 and M = 10^ (thick lines), and HFT theory 



a given — 1/\/2j/, see eq. ([4|, decreases as grows 
so larger system sizes are needed to observe convergence 
to the macroscopic limit. In addition, finite-size effects 
related to the number of clones M used for the sam- 
pling become an issue in this limit [5]. In any case, the 
sharpening of G{d) as i' increases shows that large dissi- 
pation fluctuations are strongly supressed in this regime. 
In the strongly-inelastic limit, i/ ^ 1, the scale — ?> 
and the system decouples effectively into two indepen- 
dent boundary shells, with energy concentrated around 
the boundary baths. This is evidenced by the optimal 
energy proflles for a given d measured for v = 10, see 
bottom panel in Fig. [3] in contrast to the behavior ob- 
served for <C 1, see Fig. [2] Moreover, a simple scaling 
argument shows that G{d) \fvT(dj \fv), with T(^z) 
some scaling function ,13j . 

Recently a similar hydrodynamic theory has been de- 
veloped to study large fluctuations in driven dissipative 
media |17j . but its predictions do not compare well with 
numerical results of the dissipative KMP model [13]. 
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